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Charge correlations in lattice gauge calculations suggest that up, down and strange charges move 
independently in the QGP (quark-gluon plasma), and that the density of such charges is similar 
to what is expected from simple thermal arguments. Here, we show how specific elements of the 
charge-charge correlation matrix in the QGP survive hadronization and become manifest in final- 
state charge-charge correlation measurements. 
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Identifying partonic degrees of freedom in the QGP 
(quark-gluon plasma) is a principal goal for the rela- 
tivistic heavy ion programs at RHIC (Relativistic Heavy 
Ion Collider) and at the LHC (Large Hadron Collider). 
From the strong collective flow observed at RHIC, it 
is clear that the QGP interacts strongly and behaves 
like a liquid [1]. The low viscosity suggests that col- 
lisional widths of partons are similar in magnitude to 
their energy, so quarks and gluons are not particularly 
good quasi-particles. Nonetheless, they can still serve as 
the fundamental carriers of charge. Lattice gauge cal- 
culations confirm the picture that up, down and strange 
charges are transported as single charges, in contrast to 
the hadronic world where they are transported in either 
charge-anticharge pairs (mesons) or in groups of three 
charges (baryons) [5J [3] . Insight from lattice calculations 
comes from charge-charge correlations. 



Xab = {QaQb)/V, 



(1) 



where within the volume V, Qa is the net number of 
up, down or strange quarks. Here, we consider systems 
with zero net charge, otherwise one would subtract terms 
{Qa){Qb)- If each charge carrier had one unit charge and 
moved independently, as in a parton gas, the off-diagonal 
elements would be zero and the correlation would be 



(parton gas) 



Sab (na 



(2) 



For example, n„ is the density of up quarks and n-a is 
the density of anti-up quarks. The non-zero correlation 
comes from charges being correlated with themselves. 
This contrasts with a hadronic gas, 
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where qa.a is the charge of type a on the hadronic species 
a. Mesons tend to make the off-diagonal elements nega- 
tive while baryons have the opposite effect. Within the 
last few years lattice calculations have provided Xab as a 
function of temperature, and shown the transition from 
a hadron gas to a QGP [2 |3]. For temperatures above 
the critical region, the off-diagonal elements largely dis- 
appear and the magnitude of the diagonal elements ap- 
proach 90% of the expected value for a parton gas. In 



this letter we show how Xab from the QGP phase becomes 
manifest in experimental measurements. 

Local charge conservation precludes correlations from 
lattice calculations being directly compared to those in 
experiment. Strictly speaking, the net charge in the com- 
plete volume is fixed and does not fiuctuate, whereas in 
lattice calculations one assumes interaction with a par- 
ticle sink, i.e., the grand canonical ensemble. However, 
charge correlations binned as a function of relative coor- 
dinate can be analyzed, and depend on Xab as calculated 
in a grand canonical ensemble. For the partonic stage in 
a heavy-ion collision, we assume the following form for 
the binned correlation. 



9abiv,'r) = {Pa{0)pbiv)), 
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Since quarks are relatively light (m < T), the number 
of quarks is roughly fixed between the initial thermaliza- 
tion of the QGP and hadronization. Balancing charges 
then diffuse away from one another with a separation 
described by ctqgp. The delta function represents the 
correlation of a charge with itself, and if the balancing 
charges were allowed to spread over an infinite volume 
it would be the sole term and one would reproduce the 
parton gas value for Xab in Eq. ([2|. The relative coor- 
dinate along the longitudinal axis is the spatial rapidity, 
rj — tanh~ (z/t), and is known as the spatial rapidity, 
and is chosen instead of the coordinate z to better ac- 
count for the longitudinal expansion, which for a boost 
invariant system is Vz = z/t. In a boost invariant system 
correlations depend only on the relative spatial rapidity 
rj. The width cqqpIt) is then frozen in the limit of zero 
diffusion. From this point forward, charge and number 
densities are per unit 77, and from global charge conser- 
vation, J dr] gabiv) — 0- Fo^ an adiabatic expansion the 
entropy and quark numbers per unit 77 are approximately 
fixed, and in the gas limit Xab would stay constant during 
the evolution of the QGP. The proper time r = y/t'^ — z^ 
is the time measured by an observer moving with the fluid 
from the origin z = t = 0. 

More generally, the strength of the Gaussian in the ex- 



pression for Qabiv) in Eq. (U| is xib ; which allows one 
to compare to lattice calculations. Unfortunately, mea- 
surements are not made before hadronization, but only 
afterward, when particles are on their final-state trajecto- 
ries. To connect to experiment one must overcome three 
challenges: 

a) A second wave of quark-antiquark pairs are created 
or destroyed at hadronization. 

b) Hadrons are measured ~ not quarks. 

c) One can measure only momenta, not spatial ra- 
pidity. Furthermore, strangeness violating decays 
affect the conservation rules. 

To overcome (a) one uses the fact that gabiv) cannot 
change suddenly except at r; '^ due to local charge con- 
servation. Additionally, the correlation afterward must 
integrate to zero. Assuming hadronization is relatively 
sudden, the correlation after hadronization should have 
the form. 
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The prime on (^'(had) denotes that the correlation of a 
particle with itself is subtracted out. The width chad 
represents the separation of those balancing charges cre- 
ated or destroyed at, or just after, hadronization and 
should be significantly smaller than ctqgp. The quantity 
Xab is defined in Eq. ([3|, and since the yields Ua are 
measured, Xab i^ ^ known quantity. The form for g'^^^ 
in Eq. (l5| comes from freezing the long-range part, while 
constraining gab to integrate to zero, or equivalently, en- 
forcing g'^^ to integrate to -xi&^°^- Thus, g'^^{r]) de- 
pends on: Xab ' "which can be taken from Eq. i2h 
or even directly from lattice calculations, Xab ' which 
can be extracted from measured hadronic yields, and the 
width Ctqgp- 

For (b) it was shown in [4] how to use g'(^^^) to gen- 
erate correlations between hadronic species, 

Go^piv) = ([rv(0) - na(0)][n^(77) - n^(r,)]), (6) 

where the Greek subscripts a and /3 refer to spe- 
cific hadronic species. To derive Gap{ri) from g'abiv) 
one assumes that correlations are distributed thermally 
amongst species, i.e.. 



{naiO)np{ri)) = (n„(0)) (n^(?7))e 



tJ.ab(v)'lc.aqi3.b 
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where ^ab is a Lagrange multiplier responsible for con- 
straining g^^. Because g'^^^ can be determined from G^p, 



a'abiv) = (l/4)X!^"/3('?)'?«,ag^,fc 



(8) 



one can write g' in terms of G, then use Eq. (|6| to write 
G in terms of correlations of densities, which using Eq. 
([t]) above to ultimately express g'{r]) in terms of //(ry) 
and measured densities (na). In [1] it was shown how to 
invert the expression and find fi in terms of g' assuming 
/i is small. This led to an expression for G in terms of g' , 



Gapijl) =4^(nc,(0))<7„,„ 
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Inserting the exprssion for g' in Ec^. ^ into Eq. (|9| one 
sees that G^p [r]) between any two hadronic species is de- 
termined by Xab ' ^he measured yields of hadrons, and 
two unknown parameters ctqgp and cthad- One doesn't 
expect significant sensitivity to cthad since it should be 
smaller than the thermal spread, so effectively the un- 
knowns in calculating Gapir]) are Xat ^n'^ ctqgp. 

Finally, to overcome (c), the correlations in spatial 
rapidity must be mapped to those in regular rapidity 
y = tanh^ {"^z)- Thermal smearing in the mapping of 
Gapiv) ~^ Gap{y) and the decays of strange hadrons 
are taken into account with a blast- wave model [i[ . The 
blast-wave model was based on a gaussian distribution 
of transverse flow velocities (with a variance of 0.6c) and 
a fixed kinetic breakup temperature (120 MeV). Param- 
eters were set to match mean transverse momenta for 
pions and protons measured by the STAR collaboration 
at RHIC JL2J. After generating correlated pairs with the 
weights described by Gap, the decays of hyperons and 
neutral kaons were simulated. The final products along 
with their weights were binned to generate Gap{y)- De- 
tails of the method can be found in [3] , though the blast- 
wave model used in this study is slightly different since it 
uses a Gaussian profile. The fit to baryon yields was ac- 
complished by assuming a chemically equilibrated sam- 
ple with a temperature of 165 MeV, combined with a 
suppression factor for baryons of two thirds to roughly 
match proton yields measured by PHENIX [T^] . The de- 
fault calculation assumes that the density of all quarks 
during the QGP is 0.85 times the final-state density of 
hadrons, and that the density of strange quarks is 92% 
that of up or down quarks. This last assumption cffec- 
tively sets the default Xa^ ■ Aside from a factor of the 
yields, Gap{y) is equivalent to charge-balance functions 
which have been studied in depth |3{S] and measured for 
a few species [5HTT]. 

The goal of this letter is to illustrate how Gap{y) is sen- 
sitive to specific elements of Xa^ and ctqgp • Figure ll 
shows how the width of the correlation 0^+^- depends 
on CT(QGP). For pions, roughly two thirds of the corre- 
lation comes from the hadronization component (with a 
characteristic spread of cthad, while approximately one 
third comes from the QGP component. This is expected 
given that well over half of the up and down quarks 
should be created at hadronization [5]. The fraction can 
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FIG. 1. The correlation for tt+tt" scaled by multiplicity and 
separated by the QGP and hadronization contributions in the 
upper panel. The tail of the correlation function is dominated 
by the QGP portion, and is more pronounced for larger (jqgp 
as is seen in the lower panel. 
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be understood through entropy arguments, which suggest 
that the number of quarks in the QGP and final-state 
hadrons are roughly equivalent, but since hadrons have 
multiple quarks, the number must more than double at 
hadronization. In Fig. [TJcthad was set to 0.2, while re- 
sults for three different values of cr(QGP) are shown, the 
default value of 1.0, along with 0.5 and 1.5. Currently, 
there are no good experimental handles on ctqgp, so an 
analysis like this is required. One can make estimates 
of the diffusive width with crude diffusion models [5]. 
Such estimates are > 0.5, but may understate the case 
since the balancing charges created at the formation of 
the QGP may also be separated due to the tunneling as- 
sociated with the breaking of color flux tubes. The tail 
of the correlation in Fig. IT] is clearly sensitive to ctqgPj 
but will be difficult to distinguish unless one has a de- 
tector with a wide acceptance in rapidity. The STAR 
detector at RHIC and the ALICE detector at the LHC 
cover ranges of relative rapidity close to 2 units, which 
may make the determination of itqgp difficult. Although 
the CMS and ATLAS detectors do not have particle ID 
for a large rapidity range, they can identify charges, and 
cover well over three units of rapidity. This should be 
sufficient for determining ctqgp at LHC energies. 

To investigate the sensitivity of the correlations to spe- 
cific elements of Xai i several correlations are plot- 
ted in Fig.s [2J [3] and |4] If the diagonal elements are 
reduced by a factor of two, the dip in the pp correla- 
tion disappears as seen in the middle panel of Fig. [2] 
This would be the case if the QGP had half the ex- 
pected number of charge carriers. If the ss component of 



FIG. 2. Correlations involving other hadronic species are 
more effective at separating the QGP and hadronization com- 
ponents because the two components might contribute with 
opposite signs. Such is the case for the pp and pK~ correla- 
tions shown in the upper panels of Figs. [2] and [3] respectively. 
The dip in the pp correlation and the switch of signs in the 
pK~ correlation only exist if there are competing correlations 
of opposite sign. Proton-antiproton correlations test the idea 
that charge production comes in two separate waves, because 
the hadronization and QGP contributions have opposite signs 
as seen in the upper panel. The middle panel illustrates the 
sensitivity to reducing the diagonal elements of Xat ^y ^ 
factor of two, and the effect of adding large off-diagonal com- 
ponents as in Eq. (fToll is displayed in the lower panel. 
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were halved while leaving the uu and dd compo- 
nents unchanged, the K'^K" correlation, shown in Fig. 
|4J becomes much narrower. This is expected, because 
in default case the hadronization contribution is small, 
and the correlation is dominated by the longer-range 
QGP contribution. This derives from the net number 
of strange quarks in the QGP being close to what one 
expects from an equilibrated hadron gas. If the QGP 
contribution is reduced, the hadronization contribution 
grows and results in a much narrower correlation. 

The off-diagonal elements of Xa^ also affect final- 
state correlations. As an extreme case, we consider a 
case where the QGP is dominated by correlated quark- 
antiquark pairs. In that case there are large off-diagonal 
elements, and if the quarks are equally well paired with 
quarks of all species, one would expect (assuming isospin 
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FIG. 3. The hadronization and QGP components of the 
pK~ correlations, shown in the upper panel, have opposite 
signs, which results in a correlation that changes sign. This 
would represent dramatic evidence of the two- wave nature of 
charge production. The pK~ correlation is sensitive to the 
off-diagonal elements of x i a-s can be seen by comparing 

the upper panel where there are no such elements to t he lower 
panel which uses the elements described by Eq. ( 10 1. 



symmetry between up and down): 
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(10) 



energy scan that covers the critical region for QCD. 
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FIG. 4. Since few strange quarks are produced during 

hadronization, the hadronization component for K'^ K~ cor- 
relations is small as shown in the upper panel. If the QGP 
had half the strangeness content, the hadronization compo- 
nent would be larger and one would have a much narrower 
correlation as seen in the lower panel. 



For this case of strong off-diagonal elements, both pp cor- 
relations (lower panel of Fig. p| and pK~ correlations 
(lower panel of Fig. Isl) are substantially affected. The ef- 
fect comes from the fact that observing a charge does not 
imply observing a unit of baryon charge since all quarks 
are matched by antiquarks. Thus, one doesn't expect to 
find a balancing charge far away and the QGP compo- 
nents to the correlations featuring at least one baryon 
nearly vanish. 

The correlations Gap {y) provide the means to connect 
Xabi measured on the lattice, to observables in heavy ion 
collisions. The four unknown elements of Xab (using 
isospin symmetry) along with the unknown width, ctqgPj 
can be extracted from analysis of Gap , which for the set 
of protons, kaons and pions, encompass six independent 
distributions, some with non-trivial shapes. The simple 
two-wave picture of charge production used here is only 
approximate and assumes that Xab stays constant during 
the QGP, but the ideas can be incorporated into a model 
with continuous evolution of x- Although the calcula- 
tions presented here can only be applied at high ener- 
gies, where there are roughly equal numbers of particles 
and antiparticles, the methods can be extended to in- 
corporate non-zero net charge. The methods could then 
address the evolution of the correlations through a beam 
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